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ABSTRACT 



( Dra (HIP 75458) is a well-known example for a K giant hosting a substellar companion since its discovery bv lFrink et alj ( 120021) . 
We present radial velocity measurements of this star from observations taken with three different instruments spanning nearly 8 years. 
They show more clearly that the RV period is long-lived and coherent thus supporting the companion hypothesis. The longer time 
baseline now allows for a more accurate determination of the orbit with a revised period of P = 5 1 1 d and an additional small linear 
trend, indicative of another companion in a wide orbit. Moreover we show that the star exhibits low amplitude, solar like oscillations 
with frequencies around 3-4 d"' (34.7-46.3 yuHz). 
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1. Introduction 

Up to now long period radial velocity (RV) variations have been 
detected in severa l K giants, e.g. P Gem (iHatzes et al.l l2b06t 
iReffert et al. 2006), HD 47536 (Setiawan et al."2003"), HD 13189 
(Hatzes et al. 2005J and t Dra ( Frink et a l. 2002). Rotational 
modulation can be excluded as a cause of these variations for 
most of these giants due to lack of variations in photome try or 
line bisectors with the RV period jHatzes & Cochranll998l) . The 
most likely interpretation is orbiting, substellar companions. 

Since the progenitor stars to planet hosting giant stars can 
have masses significantly larger than IMq, these giant stars off'er 
a way to investigate planet formation around massive stars. In 
their evolved phase the Doppler induced eff'ects of a planet are 
easier to detect than in their main-sequence phase, due to cooler 
effective temperatures and smaller rotational velocities. 

But for evolved stars the determination of the mass is more 
difficult. The most practical way for determining the stellar mass 
is the use of evolutionary tracks, but those tracks converge in the 
region of the giant branch for stars over a wide range of masses 
and thus makes the derived masses uncertain. 

Another possibility to measure the stellar mass is via astero- 
seismology. This requires the investigation of stellar oscillations. 
Photometric and RV variations with short periods from hours to 
days ha ve already been de t ected in some K and late G giants (e.g. 
a Bo o: iRetter et al." 2003; Tarra nt et al.ll2007l: a Ari : 'Kim et all ^ (jESA'l997) 
2006t e Qph: de Ridder et al. 2006t f H y a: iFrandsen et a l. 20021 " |van Leeuw^[2007l) 
Stelloet all 120061; /3 UMF lTarrant et all l2008l) and are consis- dRichichi et al j2005h 



nation with interferometric measurements of the angular diam- 
eter, an estimation for the stellar mass completely independent 
from evolutionary tracks. In a similar way we will do this here 
for t Dra. 



2. Stellar Parameters 

In Table [T] we summarize some direct measurements for the 
K2III star t Dra. The impro ved parallax from the Hipparcos cata- 
log bv lvan LeeuwenI ( l2007h is given. Angular diameter estimates 
based o n spectrophoto metry are available from the CHARM2 
catalog (iRichichietal.ll2005i) . 



Table 1. Stellar parameters of t Dra. 



Parameter 



Value 



Spectral type" 
V" [mag] 
Mv [mag] 
L[Lo] 

Parallax'' p [mas] 
Distance d [pc] 



K2III 
3.29 
0.81 

64.2 ±2.1 
32.23 + 0.1 

31.03 ±0.1 



Angular Diameter'^ 6 [mas] 3.73 ± 0.04 



Radius R [Rq] 



12.38 ±0.17 



tent with solar-like p-mode oscillations. Such oscillations have 
also been measured for ensembles of red giants with photometry 
of st ar rich regions or clusters (e.g. Gilliland ' 20081 IStello et al.l 
I2OO811) which is a very efficient way to perform such asteroseis- 
mologic studies. 

Solar like oscillations have also been discovered in planet 
hostin g main seq uence stars (e.g. jj. Arae: JBazot et al. 20051 
I Hor: IVauclair et al. 2008) and in the planet hosting K giant 
jS Gem for which [Hatzes & Zechmeisteii (l2007l) gave, in combi- 



The quantity which cannot be measured for single stars di- 
rectly is the mass, which always relies on evolutionary tra cks. 
With the online to oQ from Girardi (see Ida Silva et a n i2006l and 
iGirardi et al.l 12000 for a description), we derived some values 
for M, \ogg and R. This tool uses as input parameters the visual 



http://stev.oapd.inaf.it/cgi-bin/param 
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Table 2. Stellar parameters of t Dra from literature (T, log g and [Fe/H]) and derived with evolutionary tracks (M, log g and R). 



Reference 


r[K] 








[FeAl] 


M[Mo] 


Ic 


ig^ 




RlRe] 


Soubiran et al. (2008) 


4552 




2.96 




0.16 


1.40 ±0.24 


2.40 


±0. 


10 


11.82 ±0.63 


Prugniel et al. (2007) 


4543 




2.88 




0.19 


1.41 ±0.23 


2.40 


±0. 


10 


11.94 ±0.62 


Hekker & Melendez (2007) 


4605 




2.96 




0.07 


1.39 ±0.24 


2.45 


±0. 


10 


11.19 ±0.58 


Santos et al. (2004) 


4775 ± 


113 


3.09 ± 


0.40 


0.13 ±0.14 


1.71 ±0.38 


2.61 


±0. 


12 


10.34 ± 0.40 


Gray et al. (2003) 


4526 




2.64 




0.11 


1.31 ±0.24 


2.37 


±0. 


11 


11.97 ±0.65 


Prugniel & Soubiran (200V) 


4491 




2.57 




0.06 


1.24 ±0.24 


2.32 


±0. 


10 


12.23 ± 0.65 


Cenarro et al. (2001, 2007) 


4498 




2.38 




0.05 


1.24 ±0.24 


2.33 


±0. 


10 


12.15 ±0.66 


Allende Prieto & Lambert (1999) 


4466 ± 


100 


2.24 ± 


0.35 




1.05 ±0.36 










Mc William (1990) 


4490 




2.74 




0.03 ±0.11 


1.23 ±0.24 


2.32 


±0. 


10 


12.22 ± 0.66 


Williams (1974) 


4530 ± 


100 


2.60 ± 


0.25 


0.29 ± 0.20 


1.40 ±0.23 


2.39 


±0. 


10 


12.06 ± 0.64 



magnitude V, parallax p, temperature T and metallicity Fe/H. 
For the latter both quantities, T and [ Fe/H], we adopt the values 
in Ta blelHas given by several authors. lAllende Prieto & Lambert! 
(Il999i) gave no metallicity but derived themselves a mass of 
M = I.OSMq which was u s ed for the calculation of the com- 
panion mass in iFrink et aP (l2002h . This value is at the lower 
mass range resulting from Table|2](1.05 - I.7IM0). This shows 
how much the tracks depend on temperature and metallicity. 
Furthermore, Table |2] allows a cross check of the log ^-values 
derived from spectroscopy and from evolutionary tracks. The 
log ^-values derived from spectroscopy are always higher, which 
illustrates the difficulty of reliable mass determinations for K gi- 
ants. 



3. Observations 

For our analysis we have RV measurements for i Dra from 
three independent instruments (Table [3] and Figure [1]). The data 
set from the 0.6m CAT (Coude Auxiliary Telescope) with the 
Hamilton Echelle Spectrograph at Lick Observatory provides 
nearly 8 years and thus the longest time baseline for the orbit 
determi nation. The data up until March 2002 were already pub- 
lished in lFrink et al.] (l2002h . 

At the Thuringia State Observatory (TLS - Thiiringer 
Landessternwarte) l Dra was observed on 43 nights spanning 
200 d. The Alfred-Jensch 2m telescope is equipped with a Coude 
Echelle Spectrograph with a wavelength coverage from 4630 
to 7370 A. These observations primarily were carried out to 
look for stellar oscillations with short periods. For this pro- 
gram on some nights up to 20 spectra (even 6 1 on the night of 
JD=2454193) were taken (Figure|3]l. 

With a rather short time baseline of 8 days some spectra were 
obtained in June 2005 with the 2.7 m telescope at the McDonald 
Observatory and the "2dcoude" spectrograph. This instrument 
provides a wavelength range from 3600 A t o 1 yum. The RV s 
were calculated using the "Austral" program (lEndl et al.ll2()00l) . 
The whole data set is plotted in Figure]?] 

All three instruments have a resolving power around R - 
^ = 60 000 and utilize an iodine absorption cell for the wave- 
length reference. Table ]3] lists the time coverage, time span T, 
number of spectra and average precision of the radial veloc- 
ity measurements cr^v These measurement errors are internal 
errors due to instrumental effects recorded with the iodine lines. 
Comparing with point to point scatter, e.g. in Figure]3] the errors 
appear reasonable. 



Table 3. Journal of observations. 



Data Set 


Coverage 


T[d] 


N 


o-Rv [m/s] 


CAT 


2000.39-2008.38 


2918 


147 


4.3 


TLS 


2006.99-2007.59 


221 


280 


3.3 


McD 


2005.44-2005.46 


8 


62 


3.9 



4. Orbital Solution 

The parameters of the orbit were determinated by weighte d least 
squares (x^-Fit) with the differential correction method CSternd 
11941 ). We fitted simultaneously five Keplerian orbital elements, 
a linear trend and three offsets for the combined data sets as the 
measurements give relative velocities and the three instruments 
have different zero points. When weighting the measurements 
one has to take into account that not only the measurement errors 
introduce an error into the parameter determination but also the 
jitter due to pulsations: = cr^^ . + crj,. The stellar oscillations 
are discussed in the next chapter but we mention here that crp 
is of the order of 10 m/s. We adopt this value which is larger 
than the typical measurement error and therefore leads to a more 
equal weighting of all measurements in our fit. 

The resulting orbital elements are listed in Table ]4] and the 
RV curve is plotted in FigurelT] In the old solution the short time 
baseline and the unrecognizable line ar trend had resu lted in an 
overestimated period (P - 536 + 6 d: lFrink et al.ll2002h . The lin- 
ear trend is-13. 8 + 1.1 m/s/yr and may be part of a longer period 
caused by a further companion. For the calculation of the com- 
panion mass_weassumed a stellar mass of I.4M0 (33% higher 
than in lFrink et all2002h . 

Table 4. Orbital parameters for the companion to l Dra. 



Parameter 


Value 


Period P [d] 


510.88 ± 0.15 


Amplitude K [m/s] 


299.9 ± 4.3 


Periastron time Tq [JD] 


2452013.94 ± 0.48 


Longitude of periastron ai [°] 


88.7 ± 1.4 


Eccentricity e 


0.7261 ± 0.0061 


Linear trend [m/s/yr] 


-13.8 ± 1.1 


Mass function" /(ra) [Mq] 


(4.64 ± 0.33) • 10-' 


Semi-major axis a [AU] 


1.34 


companion mass m sin i [Mj^p] 


10.3 



" /(m) = = JL,(Kyrr^f 



Figure ]2] illustrates the radial velocity phased to the orbital 
period. The total remaining scatter is 13.9 m/s for all data sets 
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Fig. 1. Radial velocity measurements for l Dra from three data sets and the orbital solution as solid line. The lower panel shows the 
residuals after subtraction of the orbital solution and the linear trend. (See onhne version for a colored figure.) 



(and for the individual sets, which cover very different time 
scales: 15.4m/s (CAT), 13.9m/s (TLS) and 11.4m/s (McD)). 
This scatter is a factor of 3-4 greater than the measurement er- 
rors given in Table [3] and is largely due to stellar oscillations as 
we will explain below. 

5. Short Period Oscillations 

For the analysis of variations on short time scales we subtracted 
the trend and orbit from the data sets. The remaining residuals 
during some nights for the TLS and McD data sets are shown in 
Figures |3] and m respectively. Clearly one can see RV variations 
with an amplitude of few lOm/s. 

We performed the frequency analysis of the orbit residuals 
for each data set individually due to their different time sam- 
pling. For this we used t he prewhiteni ng procedure provided by 
the program Period04 ^Lenz & Bregeiil2004l) . A Fourier trans- 
formation was performed for the data, the dominant frequency 
extracted and subtracted from the data. This procedure can be 
repeated on the subsequent residuals until no further significant 
frequencies can be found. The prewhitening has the advantage 
that confusing alias peaks (here mainly 1 day aliases) of the ex- 
tracted frequency are removed from periodograms. 



Figures |5] and |6] show the periodograms in each prewhiten- 
ing step for the TLS and McD residuals, respectively. The TLS 
residuals show peaks at 0.040 d ' and 0.004 d"', which seem 
to be 1 month-aliases (A/ - 0.036 d ') to each other. Neither 
frequency is found in the other data sets. However, the corre- 
sponding long periods (P = 25 d and P - 250 d) are not of 
further interest for our investigation of short period oscillations. 
Subtraction of the lower frequency (/o = 0.004 d"', 12.33 m/s) 
flattens the TLS orbit residuals and removes the 1 day aliases 
of /o (second panel in FigurelSj. This prewhitening step reveals 
effectively an excess at frequencies around 3.8 d '. Extraction 
of the two frequencies /i=3.45d"' (4.77 m/s) and /2=4.23d"' 
(4.26 m/s) lowers the excess power considerably. A third fre- 
quency /3 =3.75 d"' (3.86 m/s) is already near an amplitude sig- 
nal to noise ratio of S/N=4 (the mean noise level is 0.865 m/s 
in the range o f 10-2 d~'). This threshold, as suggested by 
Kuschni g et al.l ( Il997l) . can be used as a criterion for stopping 
the prewhitening procedure. 

The fit with the frequencies is also drawn in Figure[3] Nights 
with many data points fit well, while some discrepancies exist for 
nights with sparse measurements. This may indicate that more 
frequencies are present or that the modes have a finite lifetim e 
due to stochastic excitation mechanisms (e.g lBarban et al.l2007h . 
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Fig. 2. Radial velocity measurements for l Dra from three data 
sets phased to the orbital period. The lower panel shows the orbit 
residuals. (See online version for a colored figure.) 
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Fig. 3. Radial velocity during several good nights measured at 
TLS (orbit subtracted). The solid curve is the fit with the TLS 
frequencies (/o, /i and /a). 



For the McD orbit residuals an excess of power (upper panel 
in Figure |6l) can be seen in the same region as for the TLS data. 
The peaks at frequencies >15 d ' in the McD data set are due to 
the sparser sampling leading to a stronger aliasing at higher fre- 
quencies. Extracted are two frequencies, /i =3.81 d"' (8.45 m/s) 
and /2=4.03 d"' (8.90 m/s), which are also illustrated along with 
the periodogram of the remaining residuals in the third panel of 
figure |6] 
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Fig. 4. Radial velocity during 8 consecutive nights measured at 
McD (orbit subtracted). 



The higher amplitudes of the frequencies in the McD data 
may be an effect of the short time baseline and the finite mode 
lifetime of solar-like oscillations seen in main seq uence stars and 
red giants (IStello et al.ll2006l; [Carrier et al.ll2007l) . 
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Fig. 5. Amplitude periodograms of the TLS RV data at each 
prewhitening step. The fourth panel shows the periodogram of 
the residuals in comparison with the extracted frequencies and 
the last panel shows the window function. The inset is an en- 
largement where the one month alias can be seen. 



The CAT residuals show a 555 d period, which is close to 
the orbital period. However the data set has insuflicient time res- 
olution to look for short period variability (< 1 day) since its 
Nyquist frequency is only 0.5 d"'. 
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A more accurate mass estimation could be done with the 
scaling relation for the frequency splitting: 



7 8 9 10 11 12 13 14 15 16 17 18 19 20 
Frequency [1/d] 

Fig. 6. Amplitude periodograms of the McD RV data at each 
prewhitening step. The third panel shows also the extracted fre- 
quencies and the last panel shows the window function. 

6. Discussion 

We compare our r esults with the scaling relati ons for solar-like 
oscillations from iKjeldsen & Bedding! (Il995h . The oscillation 
velocity amplitude Vqsc and the frequency /max of the strongest 
mode of the 5-minute-oscillations in the sun scales to other stars 
as: 



^osc — 



M/Me 



0.234 m/s 



/m 



62 d" 



I {TISlllK) 

Vosc/ms- 



3.5 



(1) 



(2) 



These relation s are valid for a wide range of convective stars 
jBed ding & Kieldsen 2003) and seem to be valid for other K gi- 
ants (iStello et al..,2008.) . For a power excess frequency of around 
3.8 d"', which is present in two data sets, the expected ampli- 
tude should be 6.75 m/s according eq.(|2| (based onT - 4490 K, 
[McWilliam (1990)). This agrees with the excess power seen 
in the periodograms and is comparable to the amplitude of fx 
(4.77 m/s (TLS) and 8.90 m/s (McD)). So this frequency satisfies 
the relation between amplitude and frequency typical for solar- 
like oscillations. 

Under the assumption that the simple scaling relations are 
valid, the combination of both equations to LjM oc T^-^ /f^^^ im- 
plies a luminosity to mass ratio of L/M - 29Lq/Mq for this fre- 
quency. Using L - 64.2 + 2.IL0 from the Hipparcos catalog one 
would estimate a mass of 2.2Mq (or with /^ax °^ gl Vt a \ogg 
of 2.54). Even if we roughly estimate an uncertainty of ±O.6M0 
due to strong 1 day aliasing (i.e. derived from A/|nax = ±1 d"'), 
this value seems too high in comparison with Table|2] For a mass 
of I.4M0 the expected frequency is 2.4 d '. We cannot exclude 
that we miss such a mode as 0.42 d periods are difficult to detect 
in short nights with single site observations. 



A/,) = 



M/Mo 



11.66 d 



(3) 



since the radius can be obtained from angular diameter measure- 
ments. With the assumed mass of I.4M0 for l Dra the expected 
frequency splitting is around 0.33 d"^ But the examination of 
the frequency splitting needs more effort and much more data 
for the identification of many modes. Unfortunately, our current 
data sets are not suitable for this kind of analysis. 



7. Conclusion 

We revised the orbit solution for the companion of t Dra. The 
orbital period is P = 5 1 1 d, some what lower than the value in 
the discovery paper bv Frink et all (1200 2). Furthermore there is a 
linear trend of -13.8 m/s/yr present, possibly caused by another 
companion. 

An excess of power around 3.8 d ' was found independently 
in two data sets, taken at different times and with very different 
sampling. The amplitude and the location of the power excess 
are consistent with solar-like oscillations. 

Our analysis of the short period oscillations indicates a 
somewhat high stellar mass of M - 2.2Mo, con siderably higher 
than the I.O5M0 of Allende P rieto & Lambert (1999) or the 
stellar masses derived from the Girardi track (1.2 - I.7M0). 
However, this is a preliminary result based on limited data. 
Assuming 1 .4M0 for the mass l Dra yields a minimum mass of 
lOMjup for the companion. 

Our RV measurements for t Dra indicate that it shows multi- 
periodic stellar oscillations. This means that an asteroseismic 
analysis can yield an accurate mass. This is best done by mea- 
suring the frequency splitting in the p-mode oscillation spectrum 
which requires more measurements and better sampling than the 
data presented here. 
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